PHASE SHIFT FRINGE ANALYSIS METHOD AND APPARATUS 

USING THE SAME 



RELATED APPLICATIONS 
[0001] This application claims the priorities of Japanese Patent Application No. 
2001-23200 filed on January 31, 2001 and Japanese Patent Application No. 
2001-399179 filed on December 28, 2001, which are incorporated herein by 
reference. 

BACKGROUND OF THE INVENTION 
Field of the Invention 
[0002] The present invention relates to a phase shift fringe analysis method 
\t using Fourier transform when analyzing a fringe image by using a phase shift method, 
j-j and an apparatus using the same. In particular, the present invention relates to a 

f* phase shift fringe analysis method comprising the steps of obtaining image 

u 

in information of interference fringes and the like while shifting a phase by using a 
phase shift device such as PZT (piezoelectric device), and analyzing thus obtained 

C3 

U plurality of image data items having a fringe pattern of interference fringes and the 
?n like, thereby attaining highly accurate phase information of an object to be observed; 
' f Z and an apparatus using the same. 

Description of the Prior Art 
[0003] While light-wave interference method, for example, has conventionally 
been known as important means concerning precise measurement of object surfaces, 
there have recently been urgent needs for developing an interferometry technique 
(sub-fringe interferometry) for reading out information from a fraction of a single 
interference fringe (one fringe) or less from the necessity to measure a surface or 
wavefront aberration of 1/10 wavelength or higher. 

[0004] An example of typical techniques widely in practice as such sub-fringe 
interferometry is a phase shift fringe analysis method (also known as fringe scanning 
method or phase scanning method) disclosed in "PHASE-MEASUREMENT 
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INTERFEROMETRY TECHNIQUES," PROGRESS IN OPTICS , VOL. XXVI (1988), 
pp. 349-393. 

[0005] In the phase shift method, a phase shift device such as a PZT 
(piezoelectric) device is used for phase-shifting the relative relationship between an 
object to be observed and a reference, interference fringe data is captured each time 
a predetermined step amount shifts, so as to measure the interference fringe 
intensity of each point of the object surface, and the phase of each point of the 
surface is determined by using the results of measurement. 

[0006] In the case of a 4-step phase shift method, for example, respective 
interference fringe intensities I lf I 2 , I 3 , and I 4 at individual phase shift steps are 
expressed as follows: 
q l x (x 9 y) = I 0 (x, y){l + y{x, y )cos[^, y)J 

□ I 2 (x, y) = Z 0 (x, y){\ + r (x, y)cos[^, y)+ */2} ^ 

\1 ' 3 M= h{ x >y){ l +y( x >y) cos [<f>( x >y)+xl 

J J h {*> y) = 'o(*> yW + fa y)™it( x > y)+ V 2 ]} 

Is 

£ where 



u x and y are coordinates; 



rU 



<t> (x, y) is a phase; 



\f I 0 (x, y) is the average optical intensity at each point; and 

T (x, y) is the modulation of interference fringes. 
From these expressions, the phase 0 (x, y) can be determined as 

r iM x ^)-i2M (3) 



[0007] Though the phase shift method enables measurement with a very high 
accuracy if a predetermined step amount can shift accurately, it has been 
problematic in terms of measurement errors occurring due to errors in step amount 
and in that it is likely to be affected by disturbances during measurement since it 
necessitates a plurality of interference fringe image data items. 
[0008] As a sub-fringe interferometry technique other than the phase shift 
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method, attention has been focused on one using a Fourier transform method, for 
example, as disclosed in "Basics of Sub-fringe Interferometry," Kogaku , Vol. 13, No. 
1 (February 1984), pp. 55 to 65. 

[0009] The Fourier transform fringe analysis method is a technique in which a 

carrier frequency (caused by a relative tilt between an object surface to be observed 

and a reference surface) is introduced, whereby the phase of the object can be 

determined with a high accuracy from a single fringe image. When a carrier 

frequency is introduced without the initial phase of the object being taken into 

consideration, the interference fringe intensity i(x, y) is represented by the 

following expression (4): 

i(x, y) = a(x, y)+b(x, y)cos[2rf x x + 2^ y + <fi(x, y)] (4) 

where 



□ a(x, y) is the background of interference fringes; 

in 

U b(x, y) is the visibility of fringes; 

|pi <t> (x, y) is the phase of the object; and 



f x and f y are respective carrier frequencies in x and y directions 
expressed by: 



fy f 2. tang, f 2 ^0y 
J x : ' Jy- 



X ' I 

where A is the wavelength of light, and 6 x and 9 y are respective inclinations of the 
object surface in x and y directions. 

[0010] The above-mentioned expression (4) can be deformed as the following 
expression (5): 

i(x, y) = a(x 9 y)+c(x 9 y)exp[i(2^ + 2rf 7 )]+ c* (x, y)exp[- i(2rf x + 2rf y )] (5) 
where c(x, y) is the complex amplitude of interference fringes, and c*(x, y) is the 
complex conjugate of c(x, y). 

[0011] Here, c(x, y) is represented as the following expression (6): 

c{x ,y) J( X 'y) CX f*( X < y )\ (6) 
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[0012] When the above-mentioned expression (5) is Fourier-transformed, the 

following expression (7) is obtained: 

/(7^)=^^)+^-/ :( <-/ y ) + C-(^ + / x ^ + / y ) (7) 

where 

A( V , f ) is the Fourier transform of a(x, y); 

C( n -f x , £ -f y ) is the Fourier transform of v)exp[i(2^ + 2;/ J ; and 
C*( n +f x , £ +f y ) is the Fourier transform of c*(x 9 v)exp[_z(2^ + 2,/J. 
Subsequently, C( 77 -f x , £ -f y ) is taken out by filtering, and the peak of a 
spectrum positioned at coordinates (f x , f y ) is moved to the origin of a frequency 
coordinate system (also known as a Fourier spectrum coordinate system; see Fig. 8), 
so as to eliminate the carrier frequency. Then, inverse Fourier transform is carried 
out, so as to determine c(x, y), and the wrapped phase <j) (x, y) is determined by the 
following expression (8): 

rv ' Re[c(x,y)J 
where Im[c(x, y)] is the imaginary part of c(x, y), and Re[c(x, y)] is the real part of 



in 
»* 
\* 

m 

s 

\t c(x, y) 

fU 

Hi 



[0013] Finally, unwrapping is carried out, so as to determine the phase 3>(x, y) 
of the object to be measured. 

[0014] In the Fourier transform analysis method explained in the foregoing, the 
fringe image data modulated by the carrier frequency is subjected to Fourier 
transform as stated above. 

[0015] In general, as mentioned above, the phase shift method captures the 
brightness of an image while imparting a phase difference between object light of an 
interferometer and reference light by a phase angle which is an integral fraction 2tt 
and analyzes thus captured brightness, thereby theoretically enabling a highly 
accurate phase analysis. 

[0016] For securing highly accurate phase analysis, however, it is necessary 
that the relative relationship between a sample and a reference be displaced at a 
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high accuracy by a predetermined phase amount (very short distance). When the 
phase shift method is carried out by physically moving a reference surface or the 
like by using a PZT (piezoelectric device), for example, it is necessary that the 
amount of displacement of the PZT (piezoelectric) device be controlled highly 
accurately. However, the displacement error of the phase shift device or the tilt 
error of the object surface is hard to eliminate completely, whereby controlling the 
amount of phase shift or tilt with a high accuracy is a difficult task in practice. 
[OOlffl From such a viewpoint, the assignee has proposed a technique for 
detected the above-mentioned error resulting from the phase shift device and 
correcting the fringe image analysis according to the detection value upon carrying 
out the analysis (Japanese Patent Application No. 2000-277444). In this technique, 
the fringe image data obtained by use of the phase shift method is subjected to 



I* 

H 

|2 Fourier trans\or,ni, the carrier frequency and complex amplitude caused by 

If! 

J* fluctuati ns in w%£e\fi nt occurring between in >bj< to be ol erved nd n fei nee, 

W and the amounts di displacement and tilt of phase shift are detected according to the 

s carrier frequency a\d complex amplitude, so as to correct results determined by the 

(3 

H phase shift method\ whereby influences caused by errors in the amount of 

FU 

\n tilt/displacement of theV)hase shift amount are eliminated. 

fy [0018] The above-mentioned technique proposed by the assignee is a quite 



effective technique in that it can easily alleviate influences of errors caused by the 
phase shift device without using an expensive phase shift device. However, since 
calculations for correcting the amount of error from a predetermined step amount 
(e.g., 90 degrees in the case of 4-bucket (step) method), the time required for fringe 
analysis becomes longer. Also, even in the case where 3-bucket method is 
theoretically sufficient when no shift error exists in the phase shift device, it is 
necessary to use 4-bucket method or 5-bucket method from the viewpoint of 
accuracy, which may also hinder the time required for fringe analysis from becoming 
shorter. 

[0019] In view of the foregoing circumstances, it is an object of the present 
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invention to provide a phase shift fringe image analysis method which can eliminate 
influences caused by errors in the amount of displacement of phase shift and/or 
relative tilt amounts between the object and the reference without making the 
apparatus configuration complicated and expensive when analyzing fringe image data 
obtained by using a phase shift method, whereby the fringe analysis can be carried 
out rapidly and favorably; and an apparatus using the same. 

[0020] The present invention provides a phase shift fringe image analysis 
method comprising the steps of shifting an object to be observed and a reference 
relative to each other by using a phase shift device, obtaining fringe image data at a 
plurality of phase shift positions, and determining a phase of the object by analyzing 
thus obtained plurality of fringe image data items; 

wherein the plurality of phase shift positions are at least three phase 
positions having a given phase gap therebetween; and 

wherein positional data of the above-mentioned at least three phase 
positions are specified, and the whole or part of the fringe image data on which 
carrier fringes at these phase positions are superposed is subjected to a 
predetermined arithmetic operation so as to carry out a phase analysis and 
determine the phase of the object. 

[0021] The positional data of the above-mentioned at least three phase positions 
may be determined by a Fourier transform fringe analysis method. 
[0022] The predetermined arithmetic operation may be carried out in view of 
data concerning relative tilt between the object and the reference at the 
above-mentioned at least three phase positions. 

[0023] The data concerning relative tilt between the object and the reference 
may be determined from a difference in frequency of the carrier fringes. 
[0024] The data concerning relative tilt between the object and the reference 
may be determined from a difference in phase of the object. 

[0025] The number of phase shift positions for determining the fringe image 
data may be 3, and the phase of the object may be represented by the following 
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conditional expression (1): 



where 



where 



where 



(x, y) = arctan 



cos J 3 _(l + /?)cos£ 2 + p cos St 
sin - (1 + p)sin S 2 + P sin J, 



(1) 



cos[^(x, y)+ £ 3 ]- cos[^, y)+S 2 ] 



, and 



cos[^(x, y)+ S 2 ]- cos[^(x, y)+ § x 

L (*> y^ m )= y) + b ( x > y) c ™[ 2 *fxm x + V ym y + y)+ <j m ] 

= a(x, y) + y)cos[^(x, y ) + § m ] 

a(x, y) is the background of interference fringes; 

b(x, y) is the visibility of fringes; 

0 (x, y) is the phase of the object; and 

5 m is the phase shift amount of the phase shift device expressed by: 



f m is the phase of the phase shift device (not including the part involved 
with the tilt of the phase shift device); and 

f xm and f ym are the carrier frequencies (including the part of the error in 
inclination of the phase shift device) after the m-th phase shift expressed by: 



f = 

J xm 



2.tan<9, 



_ 2. tang,. 

•> ym - 



where 

A is the wavelength of light; 

6 xm and 0 ym are respective inclinations of the object surface upon the 
m-th phase shift in x and y directions; and 

z m is the amount of displacement of the phase shift device at the m-th 
shift position (not including the part involved with the tilt of the phase shift device). 
[0026] The phase shift fringe analysis method in accordance with the present 
invention may comprise the steps of determining a complex amplitude of a fringe 
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image by the Fourier transform fringe analysis method, and obtaining the 
above-mentioned at least three phase positions according to thus determined 
complex amplitude. 

[0027] The phase shift fringe analysis method in accordance with the present 
invention may comprise the steps of selecting a plurality of sets of at least three 
local fringe image data items corresponding to each other from fringe image data at 
the above-mentioned at least three phase shift positions, obtaining positional data of 
the above-mentioned at least three phase positions concerning each set according to 
the fringe image data of respective set, and averaging positional data of phase 
positions by a number corresponding to the number of the sets, so as to determine 
final positional data of the above-mentioned at least three phase positions. 

It [0028] The fringe image may be an interference fringe image. 

j~ 5 [0029] The present invention provides a phase shift fringe analysis apparatus 

* u for shifting an object to be observed and a reference relative to each other by using a 

l~ 

Cfj phase shift device, obtaining fringe image data at a plurality of phase shift positions, 

*3 

2 and analyzing thus obtained plurality of fringe image data items so as to determine a 

0 

!=& phase of the object; 

[U 

jy wherein the plurality of phase shift positions are at least three phase 

j[j positions having a given phase gap therebetween; and 

wherein the apparatus comprises data acquiring means for obtaining 
positional data of the at least three phase positions, and phase analysis means for 
carrying out a phase analysis by subjecting the whole or part of the fringe image data 
on which carrier fringes at these phase positions are superposed to a predetermined 
arithmetic operation. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0030] Fig. 1 is a flowchart for explaining the method in accordance with an 
embodiment of the present invention; 

[0031] Fig. 2 is a block diagram of an apparatus for carrying out the method 
shown in Fig. 1; 
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[0032] Fig. 3 is a graph for explaining characteristic features of the method in 
accordance with the present invention; 

[0033] Fig. 4 is a graph corresponding to Fig. 3 in a conventional technique; 
[0034] Figs. 5A and 5B are views specifically showing a part of Fig. 2; 
[0035] Fig. 6 is a view showing a fringe image employed for verifying the 
method shown in Fig. 1; 

[0036] Fig. 7 is a view showing the results having verified the method shown in 
Fig. 1; and 

[0037] Fig. 8 is a view showing a frequency coordinate system employed in the 
method shown in Fig. 1. 



[0038] In the following, the phase shift fringe analysis method in accordance 
with an embodiment of the present invention will be explained with reference to the 



CH [0039] This method is one comprising the steps of shifting an object to be 

s observed and a reference relative to each other by using a phase shift device, 

13 . . 

il obtaining fringe image data at a plurality of phase shift positions, and determining a 

}y phase of the object by analyzing thus obtained plurality of fringe image data items. 

in 

:f, The plurality of phase shift positions are at least three phase positions having a 
given phase gap therebetween, positional data of the above-mentioned at least three 
phase positions are specified, and the whole or part of the fringe image data on 
which carrier fringes at these phase positions are superposed is subjected to a 
predetermined arithmetic operation so as to carry out a phase analysis and 
determine the phase of the object. This embodiment relates to a case where it is 
applied to an interference fringe analysis which is a typical example of fringe image 
analysis. 

[0040] Though the phase position is not required to be a specific value in this 
phase shift fringe analysis method, it is necessary that, after the fringe image data is 
obtained, positional data of each phase position where the fringe image data was 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 



in 



drawings. 



obtained be specified. 

[0041] Namely, in the conventional phase shift fringe analysis method, 
respective fringe intensity values are plotted at predetermined phase shift positions 
(0°, 90°, 180°, 270°, and 360° here) arranged at equally-spaced intervals (at intervals 
of a = 90° here) as shown in Fig. 4, and a sine wave function f(x) passing near these 
points (a, b, c, d, e) is determined. According to such a conventional technique, the 
phase shift positions are required to be set at predetermined positions with a high 
accuracy, whereby processing for correcting them to values at these predetermined 
positions and the like are necessary in practice. 

[0042] In the phase shift fringe analysis method in accordance with this 
embodiment, by contrast, respective fringe intensity values of at least three phase 
shift positions having a given phase gap (ff„ a 2 ) are plotted as shown in Fig. 3, and 
a sine wave function f(x) passing near points (p, q, r) is obtained by using a 
computing equation which will be explained later. Here, assuming that 
f(x)=Acos(x + $)+B t 

the number of unknowns is 3, i.e., A, <j) , and B. Therefore, the function f(x) can be 
specified when positional data in at least three phase shift positions and their 
corresponding fringe intensities are known. 

[0043] An outline of this embodiment will now be explained with reference to 
the flowchart of Fig. 1, before explaining it separately in terms of specifying 
positional data at individual phase positions and calculating the phase of the object to 
be observed by using thus specified positional data. 

[0044] In this embodiment, as shown in Fig. 1, interference fringe images 
carrying form information of an object to be observed on which spatial carrier fringes 
are superposed by imparting a predetermined relative inclination between a 
reference surface and an object surface is captured by a CCD image pickup camera at 
a plurality of phase shift positions while driving a phase shift device (SI). 
Subsequently, each of thus obtained interference fringe image data is subjected to 
Fourier transform (S2), and a spatial carrier frequency (f x , f y ) is found, whereby 
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C( t) -f x , £ -f y ) is determined (S3). The spatial carrier frequency (f x , f y ) at that time 
corresponds to the total tilt formed by the predetermined inclination imparted 
beforehand between the reference surface and the object surface and the error in 
inclination generated upon driving the phase shift device. 

[0045] Subsequently, with respect to thus determined C( rj -f x , £-f y ), a 
spectrum peak located at coordinates (f x , f y ) is moved to the origin of the frequency 
coordinate system (see Fig. 8), and is subjected to inverse Fourier transform (S4) 
after eliminating the carrier frequency (including the part of the error in inclination). 
Then, a complex amplitude c(x, y) is obtained, and the relative displacement amount 
between the reference surface and the object surface is determined (S5). This 
relative displacement amount does not include the part generated by the error in 
£ inclination of the phase shift device. Namely, the relative displacement amount 

f3 represents the pure displacement amount of the phase shift device toward the 

ifl 

optical axis. 

H 

m [0046] Subsequently, without carrying out processing for moving the spectrum 

& 

5 peak located at coordinates (f x , f y ) to the origin of the frequency coordinate system, 

p 

U the above-mentioned C( r) -f x , £ -f y ) is subjected to inverse Fourier transform in a 
|y state where the carrier frequency (including the part of the error in inclination of the 
[p. phase shift device) is not eliminated, so as to determine the relative tilt amount 
between the reference surface and the object surface (S6). Thus determined 
relative tilt amount includes a predetermined inclination imparted to the phase shift 
device in order to introduce the spatial carrier frequency, and the part of the error in 
inclination of the phase shift device generated upon driving the phase shift device. 
[0047] From thus determined relative displacement amount and relative tilt 
amount between the reference surface and the object surface, each phase shift 
position (phase shift amount) of the phase shift device including the phase 
displacement amount concerning the part of the error in inclination of the phase shift 
device can be determined (S7). Then, thus determined plurality of (at least three) 
phase shift positions (phase shift amounts) are inputted into a predetermined 
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computing equation which will be explained later, so as to determine a wrapped 
phase of the object (S8). Thereafter, the unwrapped phase of the object is 
determined by a known unwrapping process (S9). 

Step of Specifying Positional Data of Each Phase Position 
[0048] This step corresponds to the processing of SI to S5 in the flowchart of 
Fig. 1. 

[0049] In general, the Fourier transform fringe analysis method can determine a 
phase from a single fringe image alone by introducing a carrier frequency (relative 
inclination between the object surface and the reference surface). When the carrier 
frequency is introduced, the interference fringe intensity is represented by the 
following expression (9): 
J ~ i{*> y) = a(x, y) + b[x, y)co^ x x + 2?f y y + j(x, y) + £ ] (9) 

j=f where 

\* a(x, y) is the background of interference fringes; 

Cfl b(x, y) is the visibility of fringes; 



*0 



0 (x, y) is the phase of the object; 



£ 

C3 

|j, £ is the initial phase of the object; and 

«y f x and f y are carrier frequencies. 

: - us. 

H [0050] As mentioned above, £ is expressed by ^ = 2 ^/ A, where A is the 
wavelength of light, and z is the displace amount of the phase shift device (not 
including the part involved with the inclination of the phase shift device). 
Therefore, the above-mentioned expression (9) can be deformed as the following 
expression (10): 

i(x,30 = fl(*,y)+c(^^ (10) 
where c(x, y) is the complex amplitude of interference fringes, and c*(x, y) is the 
complex conjugate of c(x, y): 

c(x ,y )= b ( x >y>^M x <y) + $\ . (id 

[0051] When the above-mentioned expression (10) is subjected to Fourier 
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transform, the following expression (12) can be obtained: 

/(^ff)-^ff)+c( 7 -/ ir , ff -/,) + C-^ + / J , ff + /,) (12) 

where 

A( 7? , £ ) is the Fourier transform of a(x, y); 

C( V -f x , t -f y ) is the Fourier transform of y)exp[z(2^ +2^)] ; an d 
C*( 77 +f x , £ +f y ) is the Fourier transform of c*(x,y)exp[-i(2rf x + 2^ y )]. 



[0052] As in the procedure carried out in a typical Fourier transform method, 
C( n -f x , £ -f y ) is determined by filtering. After a spectrum peak located at 
coordinates (f x , f y ) is moved to the origin of the frequency coordinate system so as to 
eliminate the carrier frequency (including the part of the error in inclination of the 
phase shift device), C( T] -f x , £ -f y ) is subjected to inverse Fourier transform (S4), 
|* whereby c(x, y) is obtained. Here, assuming that £ 0 is the initial phase of the 

^3 phase shift device, and £ t is the phase after the movement, the following expression 

IF! 

u (13): 

5 ^ )s ^ww^t (13) 

Q 

\a is obtained from the interference fringe image data at the start position. 

[y [0053] Then, from the interference fringe image data after moving the phase 

f 5 

If, shift device, the following expression (14): 

fy 

Ci{x , y) J( x >y)*MM*>y)^]} (14) 

is obtained. 

[0054] As a consequence, the following expression (15): 
is obtained. 

[0055] Therefore, the phase difference between before and after moving the 
phase shift device is represented by the following expression (16): 
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[0056] 



= tan 
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M x > y ). 



(16) 



Therefore, the displacement amount of the phase shift device is 



represented by the following expression (17): 



z = z l -z 0 



2k 



tan 



Im 








Re 


~ c ^ x >y)~ 







(17) 
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[0057] When the respective displacement amounts determined for individual 
interference fringe images obtained by a predetermined phase shift are averaged, 
displacement can be detected at a high accuracy. When detected the displacement 
of the object by the Fourier transform fringe analysis method, it is not always 
necessary to use the whole fringe image data. Displacement can also be detected 
with a sufficiently high accuracy from a part of fringe image data alone. 

Step of Calculating Object Phase 
[0058] In the following, a technique for determining the phase of the object 
according to individual phase position data of the phase shift device obtained as 
mentioned above will be explained. 

[0059] In summary, this technique is carried out such that, as shown in the 
flowchart of Fig. 1, the above-mentioned C(7?-f x , £ -f y ) is subjected to inverse 
Fourier transform in a state not eliminating the carrier frequency (including the part 
of the error in inclination of the phase shift device), so as to determine the relative 
tilt amount (including the part of the error in inclination of the phase shift device) 
between the reference surface and the object surface (S6), each phase shift position 
(phase shift amount) of the phase shift device including the phase displacement 
amount involved with the inclination of the phase shift device is determined 
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according to the relative tilt amount and the displacement amount of the phase shift 
device determined as mentioned above (S7), the wrapped phase of the object is 
obtained according to a predetermined computing equation which will be explained 
later (S8), and the unwrapped phase of the object is determined by unwrapping (S9). 
[0060] Though the technique for determining the phase in this embodiment can 
correspond to any number of buckets (3, 4, 5, 7, 11, and N buckets, etc.), only 
3-bucket, 4-bucket, and 5-bucket methods will be explained here for the convenience 
of explanation. 

[0061] In the 3-bucket phase shift method based on a given shift amount, the 
interference fringe intensity at the m-th shift position introducing the carrier 
frequency is represented by the following expression (18): 

i m (x, y,£ m )= a(x, y) + b(x, y)cos[2rf xm x + 2rf ym y + <j>(x, y) + £ m ] 



(18) 

S3 = a(x, y) + b(x, y)cos[0(x, y) + 5 m ] 

m 

? ~ where 

a(x, y) is the background of interference fringes; 



b(x, y) is the visibility of fringes; 



13 

U 0 (x, y) is the phase of the object; and 



6 m is the phase shift amount of the phase shift device expressed by: 
S m =2*f xm x + 2rf ym y + Z m (18a) 

where 

$ m is the phase of the phase shift device (not including the part involved 
with the tilt of the phase shift device); and 

f xm and f ym are the carrier frequencies (including the part of the error in 
inclination of the phase shift device) after the m-th phase shift expressed by: 



^ f - 2 " tan ^ f _ 2 ' tan ^ 

^ ' * xm ^ » J ym 



where 



A is the wavelength of light; 

9 xm and 9 ym are respective relative inclinations between the reference 
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surface and the object surface upon the m-th phase shift (including the part of the 
error in inclination of the phase shift device) in X and Y directions; and 

z m is the amount of displacement of the phase shift device at the m-th 
shift position (not including the part involved with the tilt of the phase shift device). 
[0062] The phase shift amount d m of the phase shift device in the 
above-mentioned expression (18) represents the relative phase difference between 
the reference surface and the object surface including the phase displacement 
amount involved with the pure displacement of the phase shift device toward the 
optical axis and the phase displacement amount involved with the inclination of the 
phase shift device, and is determined by the following procedure. 
[0063] First, the above-mentioned C( tj -f x , £ -f y ) is subjected to inverse Fourier 
| j transform (S6) in a state not eliminating the carrier frequency (including the part of 
^3 the error in inclination of the phase shift device), so as to determine the form of the 

in 

\* object surface in a state where the relative tilt amount (including the part of the 

m error in inclination of the phase shift device) between the reference surface and the 

& 

5 object surface is superposed thereon. Subsequently, the method of least squares is 

its. 
{ J 

used for determining a plane fitting the above-mentioned surface form, and the 
IH above-mentioned relative tilt amounts 6 xm and 0 ym are determined according to the 
If. differential coefficient of the least-square plane. 

[0064] Thereafter, the relative tilt amounts 0 xm and 0 ym are inputted into their 
corresponding expressions mentioned above in connection with expression (18), so 
as to determine the carrier frequency (f xm , f y J (including the part of the error in 
inclination of the phase shift device), whereas the displacement amount z m of the 
phase shift device determined by the procedure mentioned above is inputted into its 
corresponding expression set forth in connection with the above-mentioned 
expression (18) so as to determine the phase £ m of the phase shift device (not 
including the part involved with the inclination of the phase shift device). Further, 
thus determined carrier frequency (f xm , f ym ) and phase £ m of the phase shift device 
are inputted into the above-mentioned expression (18a), so as to determine the 
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phase shift amount 8 m of the phase shift device. 

[0065] For determining the phase 0(x, y) of the object, the following 

expression (19): 

= h ~h = cosfr(*,y) + S,]-cosy,(x,y) + S 2 ) 
i 2 _i, cos[^(x,y)+ S 2 ]-cosy>(x ) y)+ S Y ] 

will now be studied. 

[0066] From this expression, the following expression (20): 

, v cos J, - fl + p)cos + P cos J. 

v sin£ 3 -(l + p)sin£ 2 -fpsinj, 

is obtained. 

[0067] The above-mentioned procedure can determine the phase shift amount 
|f 8 m representing the relative phase difference between the reference surface and 
O the object surface. Therefore, the wrapped phase of the object form can be 

in 

\* determined by using the above-mentioned expression (20). Further, a known 

\* 

m unwrapping method can be used for determining a continuous object phase O. 

*Q 

5 [0068] In the 4-bucket phase shift method with a given shift amount, the 

a 

u following expressions (21) and (22) are used in place of the above-mentioned 



expressions (19) and (20): 



cos[0(*, y ) + S A ] - cos[^(x, y ) + S 2 ] 



p = A 2 = irv ' --w\-->'/ • ^21) 

/ 3 -i, cos[^(x,y) + ^ 3 ]-cos[^(x,y)+^] * 

, . cos£ A -cos£- + pcos£, - pCOSSi 

m(jc, y) = arctan , 2 . 3 . 1 . (22) 

v 7 sin J 4 _sin£ 2 +/?sin£ 3 -psin ^ 

[0069] In the 5-bucket phase shift method with a given shift amount, the 

following expressions (23) and (24) are used in place of the above-mentioned 

expressions (19) and (20): 

p= = 2{cos[^y) + J 4 ]-cos[^y) + J 2 ]} 

i, - 2i 3 + i 5 cosy>(x, y) + s x ] - 2 cos[^(*, y ) + £ 3 ] + cos[^(x, y) + £ 5 ] * 

. v ^ 2cos£ 4 -2cosj 2 +pcos£,-2/?cos£ 3 + /?cos£ 5 

^(x, y ) = arctan i 1 — ~ ! _ 1 — . (24) 

2sin s 4 -2sin£ 2 + psin s x -2/?sin£ 3 + psin£ 5 
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[0070] An apparatus for carrying out the above-mentioned method will now be 
explained with reference to Fig. 2. 

[0071] This apparatus is used for carrying out the method in accordance with 
the above-mentioned embodiment. As shown in Fig. 2, in a Michelson type 
interferometer 1, interference fringes are formed on the imaging surface of CCD 5 in 
an image pickup camera 4 by respective reflected luminous fluxes from an object 
surface 2 and a reference surface 3, and are inputted into a computer 7 equipped 
with a CPU and an image processing memory by way of an image input board 6. 
Thus inputted interference fringe image data is subjected to various arithmetic 
operations, and the results of operations are displayed onto a monitor screen 7A. 
The interference fringe image data outputted from the image pickup camera 4 is 
ji temporarily stored into the memory due to the processing carried out by the CPU. 
J=j The piezoelectric driving signal outputted from the computer 7 is fed into a 
* u piezoelectric driver 9, whereas a PZT (piezoelectric device) actuator 10 is driven in 
response thereto. 

£ [0072] The computer 7 subjects thus obtained interference fringe image data to 

13 

I* Fourier transform, so as to determine the (relative) phase shift amount of the phase 

\'U 

ry shift device (between the object surface and the reference surface), and an 
iTi arithmetic operation of the computing equation in which the phase shift amount is 

inputted is carried out, so as to calculate the phase of the object surface. 

[0073] Figs. 5A and 5B show two modes of the PZT (piezoelectric device ) 

actuator 10, respectively. 

[0074] As shown in Fig. 5A, the first mode comprises three piezoelectric 
elements 121, 122, 123 for supporting the backside of the reference surface 
(reference mirror) 3, whereas two lines L x , L y respectively connecting the 
piezoelectric element 121, also acting as a fulcrum member, to the other 
piezoelectric elements 122, 123 on the reference mirror having the reference 
surface 3 are arranged orthogonal to each other. A phase shift is effected when the 
three piezoelectric elements 121, 122, 123 expand/contract by the same amount. 
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When the piezoelectric element 122 expands/contracts alone, the reference surface 3 
of the reference mirror inclines in the x-axis direction so as to rotate about the y 
axis. When the piezoelectric element 123 expands/contracts alone, the reference 
surface 3 of the reference mirror inclines in the y-axis direction so as to rotate about 
the x axis. As shown in Fig. 5B, the second mode is constructed such that the 
center part of the backside of the reference surface (reference mirror) 3 is supported 
by a cylindrical piezoelectric tube 124. A phase shift is effected by an unbiased 
expansion/contraction of the piezoelectric tube 124. On the other hand, a biased 
expansion/contraction freely tilts the reference surface 3 of the reference mirror in 
x- and y-axis directions. 

[0075] For verifying the accuracy of the method in accordance with the 
J* above-mentioned embodiment, a phase analysis was carried out by the 
13 above-mentioned 3-bucket method while using three fringe images whose initial 

in 

i A phases were 0°, 76°, and 123°, respectively, in interference fringes (with 256 

I* 

Cm gradations) shown in Fig. 6 obtained when an ideal plane is observed. Thus 



a 



determined phase shift amounts were 75.9986° and 122.9973°. Fig. 7 shows the 



I* error between the results of phase analysis and the ideal plane. The maximum 
In error is about ±0.01 rad (7.9 x 10" 4 A ), which is on a negligible order. Therefore, 
|y it has been verified that highly accurate phase analysis can be carried out by the 
3-bucket method as well. 

[0076] Withdut being restricted to the above-mentioned embodiment, the 
method of the roresent invention can be modified in various manners. When 
detected the inclination of the object in the method in accordance with the 
a above-mentioned embodiment, fringe image data carrying the phase information 



from the objfeetf is acquired, the whole or part of each fringe image data is 
j subsequently ^subjected to Fourier transform, so as to determine the ^phase 
^ information of the object in the fringe image data, and the inclination of the object is 
detected according to thus determined phase information. However, for example, 
the frequency of carrier fringes in the fringe image data may be determined in place 

FK0911SEC -19- 



of the phase information in the fringe image data, and the inclination of the object 
may be detected according to thus determined frequency of carrier fringes. A 
technique for detected the inclination of fcfljLebject by determining the frequency of 
carrier fringes in fringe image data as sJor is disclosed in detail in the specification 
of commonly-assigned Japanese Patenjr Application No. 2000-397246. 
[0077] The phase shift element is not limited to the above-mentioned PZT, but 
may be those which can achieve the phase shift method by physically moving the 
reference surface or object surface or changing the optical path length by use of an 
AO element or EO element, or a transmission type element which can change the 
refractive index or the like so as to alter the optical path length when inserted into 
the reference optical path and/or observation optical path. 
[0078] Though three PZT elements are exactly positioned at the respective 
vertices of a rectangular triangle as a mode for arranging them in the 
above-mentioned embodiment, the aimed effect can be obtained if the three 
members are arranged so as to form the respective vertices of a given triangle on 
the reference mirror. 

[0079] Though the above-mentioned embodiment is explained while using 
spatial carrier frequencies as carrier frequencies, temporal carrier frequencies or 
temporal spatial carrier frequencies can also be used as the carrier frequencies of 
the present invention. 

[0080] Though the interference fringe image data is captured with a Michelson 

type interferometer in the above-mentioned embodiment, the present invention is 

similarly applicable to interference fringe image data obtained by other types of 

interferometers such as those of Fizeau type as a matter of course. 

[0081] Further, the present invention is similarly applicable to various kinds of 

fringe images such as moire" fringes and speckle fringes in addition to interference 

fringes. 

[0082] In the phase shift fringe analysis method in accordance with the present 
invention and the apparatus using the same, a phase shift device is employed for 
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phase-shifting an object to be observed and a reference relative to each other, fringe 
image data of at least three phase positions having a given gap therebetween are 
obtained, positional data of the at least three phase positions are specified when 
determining the phase of the object by analyzing thus obtained plurality of fringe 
image data items, and a phase analysis is carried out by subjecting the whole or part 
of the fringe image data on which carrier fringes at these phase positions are 
superposed to a predetermined arithmetic operation, so as to determine the phase of 
the object. 

[0083] Therefore, as compared with the prior art in which a predetermined shift 

position is fixed as a specified value, the time required for a calculation for 

correcting the amount of error from the predetermined shift position becomes 

unnecessary. Also, the accuracy in measurement can be improved by a bucket 

[2 method with a smaller number even without using expensive phase shift means, 

f* whereby the time required for fringe analysis can be shortened. 
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